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The majority of intravital studies on brain tumor in living animal so far rely on dual color 
imaging. We describe here a multiphoton imaging protocol to dynamically characterize 
the interactions between six cellular components in a living mouse. We applied this 
methodology to a clinically relevant glioblastoma multiforme (GBM) model designed in 
reporter mice with targeted cell populations labeled by fluorescent proteins of different 
colors. This model permitted us to make non-invasive longitudinal and multi-scale 
observations of cell-to-cell interactions. We provide examples of such 5D (x,y,z,t, color) 
images acquired on a daily basis from volumes of interest, covering most of the mouse 
parietal cortex at subcellular resolution. Spectral deconvolution allowed us to accurately 
separate each cell population as well as some components of the extracellular matrix. 
The technique represents a powerful tool for investigating how tumor progression is 
influenced by the interactions of tumor cells with host cells and the extracellular matrix 
micro-environment. It will be especially valuable for evaluating neuro-oncological drug 
efficacy and target specificity. The imaging protocol provided here can be easily translated 
to other mouse models of neuropathologies, and should also be of fundamental interest 
for investigations in other areas of systems biology. 

Keywords: intravital imaging, two-plioton microscopy, cellular interactions, brain tumors, multicolor transgenic 
mouse, spectral deconvolution, cranial window 



INTRODUCTION 

Over the past decade, intravital two-photon microscopy (2PM) 
(Ustione and Piston, 201 1) has become a gold standard for funda- 
mental and pre-clinical research on various pathologies including 
cancer (Perry et al., 2012), auto-immune (Niesner et al, 2013) 
and vascular (Kuijpers and Heemskerk, 2012) diseases. It has 
also been applied in the field of neuroscience to study brain 
(Trachtenberg et al, 2002), spinal cord neuronal and vascular 
networks (Misgeld and Kerschensteiner, 2006; Dray et al., 2009; 
Fenrich et al., 2012). The technique is non-invasive and imag- 
ing can be repeated in the same animal. Chronic imaging thus 
avoids the sacrifice of large cohorts of animals at different time 
points as previously required for classical histology and immuno- 
histochemistry experiments (Zhai et al, 2011; Piao et al., 2012). 
Thanks to its high spatial resolution, dynamic subcellular obser- 
vations can be performed in the brain of living transgenic mice 
(e.g., Thyl-CFP, Feng et al., 2000). Exogenous fluorescent mark- 
ers can also be used to label specific compartments such as blood 
vessels through intravenous injections of Rhodamine-dextran 
(Verant et al, 2008; Ricard et al, 2009). The outstanding imag- 
ing depth of 2PM in scattering media permits investigators to 
scan most of the subdural brain cortical layers in mice (Helmchen 
and Denk, 2005; Kobat et al., 2009). In the case of brain tumors, 
it offers unique means for studying the kinetics of tumor cell 
infiltration/migration as well as permitting access to correlated 
successive events such as vascular remodeling induced by tumor 
cell proliferation. The efficacy of anti-tumor and anti-angiogenic 



drugs can thus be assessed and their action mechanism clarified 
(Von Baumgarten et al, 201 1; Ricard et al, 2013). 

The role of the tumor microenvironment has recently become 
a hot topic in neuro-oncology (Charles et al., 2011; Hanahan 
and Weinberg, 2011). The microenvironment includes several 
populations of cells that interact directly with extracellular 
matrix proteins. All these components should ideally be moni- 
tored simultaneously to elucidate the key mechanisms underly- 
ing tumor progression. However, current imaging protocols are 
designed to study at most 3 parameters in the same intravital 
imaging sessions, and are usually limited to the study of tumor 
cells and vasculature only (Winkler et al., 2009; Von Baumgarten 
et al., 201 1). Glioblastoma multiforme (GBM) is the most aggres- 
sive form of brain tumor, with poor prognosis and no curative 
treatment available to date (Deangelis, 2001; Prados et al., 2004; 
Ricard et al., 2012b). We have therefore developed an orthotopic 
GBM murine model, consisting of the graft of fluorescently- 
labeled tumor cell spheroids into the cerebral cortex of multicolor 
fluorescent transgenic mice where each cell population and com- 
partment of interest fluoresce in a specific color. We have subse- 
quently set up dedicated 2PM and image processing protocols to 
routinely observe 6 components on a daily basis in the same living 
animal without side effects. We describe the complete methodol- 
ogy and illustrate new dynamic information that can be obtained 
at micro and macroscopic scales with imaging sessions lasting 1 h 
on average. Finally, we discuss the possible implications of this 
approach for pre-clinical studies. 
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BIOLOGICAL METHODS 
ANIMAL CARE GUIDELINES 

All experimental procedures were performed in accordance 
with the French legislation and in compliance with the 
European Community Council Directive of November 24, 1986 
(86/609/EEC) for the care and use of laboratory animals. 
The research on animals was authorized by the Direction 
Departementale des Services Veterinaires des Bouches-du-Rhone 
(license D-13-055-21) and approved by the National Committee 
for Ethic in Animal Experimentation (Section N°14; project 
87-04122012). 

CELL CULTURE 

Cells from the murine GL261 GBM cell line (Newcomb and 
Zagzag, 2009), were transfected with a plasmid encoding for 
DsRed (pDsRed2-Nl, Clontech, Mountain View, USA). Cells that 
stably express DsRed were selected using Geneticin (0.5 mg/ml, 
Gibco) and cultured as monolayers in RPMI1640 -|- GlutaMAX- 
1 (Gibco 61870) supplemented with 10% Fetal Calf Serum 
(Thermo Scientific) and Geneticin (0.5mg/ml,Gibco). Cells were 
kept at 37°C in a 5% CO2 atmosphere. Confluent cells were plated 
on Petri dishes coated with 1% soft agarose to induce spheroid 
formation. 

TRANSGENIC MICE 

Connexin43-CFP (astrocytes Degen et al., 2012), Thyl-GFP 
(neurons Feng et al., 2000) and GDI Ic- YEP (reactive microglia 
and dendritic cells Lindquist et al., 2004) mice were crossed 
to obtain the triple transgenic mouse strain, Connexin43- 
CFP/Thyl-GFP/CDllc-YFP, used in our experiments. 

GLIOBLASTOMA ANIMAL MODEL 

The technology had been described in detail in Ricard et al. 
(2013). Briefly, adult mice (>7 weeks) were anaesthetized by 
an intraperitoneal injection of a mixture of Xylazine/Ketamine 
(lOmg/kg and lOOmg/kg, respectively) and positioned on a 
stereotactic frame. A 3-4 mm diameter craniotomy was per- 
formed on the left parietal bone. The dura-mater was incised by 
a 31G needle and a 250 (Am spheroid of GL261-DsRed murine 
GBM cells was then injected into the cerebral cortex approxi- 
mately 250 |xm below the brain surface. A Sephadex hemi-bead 
with a diameter that fit through the dura-mater opening was 
inserted in the injection wound and glued using histo-compatible 
acrylic glue (Cyanolit). A round glass coverslip (5 mm diame- 
ter) was then sealed on the adjacent bone and fixed to the skull 
by dental cement. Animals were allowed to recover for 15 days 
post-surgery before the first imaging session. 

ANIMAL PREPARATION FOR IMAGING 

Prior to each imaging session, mice were anaesthetized by 
an intraperitoneal injection of a mbcture of Xylazine/Ketamine 
(lOmg/kg and l.Omg/kg, respectively), injected intravenously 
with 100 |xl of a Cascade Blue conjugate dextran (70kDa) solu- 
tion (20 mg/ml in phosphate saline buffer (Sigma)) and posi- 
tioned on a stereotactic frame. Each imaging session, including 
anesthesia induction, animal positioning and acquisition lasted 
approximately 45 min-1 h. 



EXPERIMENTAL SETUP AND IMAGE PROCESSING 

We used a Zeiss LSM 7 MP two-photon microscope home- 
modified to allow animal positioning below the 20X water 
immersion objective (1.0 NA) and coupled to a femtosecond 
pulsed infrared tunable laser (Mai-Tai, Spectra Physics). After 
two-photon excitation, epifluorescence signals were collected 
and separated by dichroic mirrors and filters on 5 indepen- 
dent non-descanned detectors (NDD) (Figure 1). Images were 
acquired sequentially using an excitation wavelength tuned at 
800 nm and then at 940 nm, and were assembled into a spec- 
tral hyperstack. Gains and offsets were identical for all the 
detectors, except for the red channels whose gain was slightly 
reduced to compensate for the very strong expression of DsRed 
in tumor cells. 

Each fluorescent reporter (i = 1-5) used to selectively label a 
physiological component exhibited a footprint on all detectors 
(NDD1-800-NDD5-800, numbered #1 to #5 and NDD 1-940 - 
NDD5-940, #6 to #10). Its spectral signature on the detection 
arrays was represented by a normalized vector Si = [Qj, 
Q,io] whose coordinates Qj: (k = 1-10) corresponded to the rel- 
ative contribution of each detector k to the sum of intensities 
measured on all detectors (in percent. Figure 2). Spectral sig- 
natures were determined from regions of interest (ROIs, n = 5) 
selected manually for their unambiguous and exclusive expres- 
sion of one of the fluorescent labels of interest. Control intensities 
measured in a ROI not containing any fluorescent object were 
systematically subtracted from the signals of fluorescent ROIs in 
order to correct for background signals prior to normalization. 
For each fluorophore (i), the main contribution to the global flu- 
orescence was identified as Q^REPi. and the detector #REFi was 
thus defined as the reference detector to outline biological objects 
labeled with the fluorophore (i). Quantitative data are depicted in 
Figure 2F. 

Spectral deconvolution of the channel #REFi was finally 
obtained by subtracting the minor contributions of all other fluo- 
rophores in the channel #REFi. Corrections for each fluorophores 
j were then expressed as percentages of the fluorescence intensities 
measured on their respective reference detector #REFj. 

After determining the spectral signature S, of the fluorophore 
of interest and the signatures Sj for all other contaminating labels 
in ROIs selectively labeled with one fluorophore, 

Si = (Qj, . . . , Q refj, Q.io) with 1 < REFi < 10 and 

Sj = (q,i, . . . , q.REFj, q.io) i < REFj < lo 

spectral deconvolution of the whole image was obtained using the 
following equation for each of its pixels: 

^i.Deconv = ^i.REF " ^ '^j ' ^j.REF 

where: /i,Deconv is the spectrally unmixed image for element i. 
where: /i,REF is the reference image for element it displays the 
strongest fluorescent signal contribution for the element i. 
where: -fj^REF is the reference image for element j; it displays the 
strongest fluorescent signal contribution for the element;'. 
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FIGURE 1 I Schematic representation of the two-photon microscopy 
setup. The excitation beam is produced by a femtosecond pulsed 
infrared tunable (720-1020 nm) laser (Mai-Tai, Spectra Physics). The 
laser power is modulated by an Acousto-Optic Modulator (AOM). The 
beam is scanned in the xy direction by gaivanometric mirrors present 
in the scan head of a Zeiss LMS 7 MP two-photon microscope. The 
beam then passes through a LP690 dichroic mirror and is focused in 



the brain of the anaesthetized animal by a 20X-1.0 NA water 
immersion objective. The emitted epifluorescence is collected and 
reflected by the LP690 mirror in a non-descanned mode. The 
fluorescence is finally splitted and filtered using a set of dichroic 
mirrors and filters and collected by a set of 5 non-descanned 
detectors mounted in cascade (NDD). The characteristics of the 
dichroic mirrors and filters are depicted on the scheme. 
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FIGURE 2 I Normalized spectral signatures of each biological element on and 940 nm (NDD1-940-NDD5-940). (F) Numerical values (%) used for 

the non-descanned detectors array. (A-E) normalized (%) contribution of spectral deconvolution (red: contribution of the element of interest in its 
each element on the NDDs for excitations at 800 nm (NDD1-800-NDD5-800) reference NDD). 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



February 2014 | Volume 8 | Article 57 | 3 



Ricard and Debarbieux 



Six-color intravital two-photon imaging 



and: 

Cj,REFi 
Cj.REFj 

where: Cj^EPi is the relative contribution of element; on the ref- 
erence NDD for element i. It can be read from Sy, 
where: Cj^REPj is the relative contribution of element; in its own 
reference NDD. It can be read from Sj . 

Spectral deconvolution was performed for the 5 fluorophores 
(CFP, GFP, YFP, DsRed, cascade blue) to generate a set of five 
image stacks representative of the distribution of the following 
components: astrocytes, neurons, CDllc-positive cells, tumor 
and vasculature. A sixth element, SHG arising from dura-mater 
and peritumoral collagen, was only present in the superficial lay- 
ers of the brain and was spectrally unmixed from NDD2-940 
using the same approach (magenta pseudo color in Figure 4). All 
spectrally unmixed individual images for each biological com- 
ponent and the collagen images were finally merged together to 
obtain a corrected 6-color image (blue: vasculature, cyan: astro- 
cytes, green: neurons, yellow: CDllc-positive immune cells, red: 
tumor, magenta: collagen). Custom- written ImageJ macros were 
used to spectrally unmix the fluorescent signals. 

RESULTS AND DISCUSSIONS 

A set of five detectors was installed in the non-descanned posi- 
tion of our 2-photon upright microscope. Filter cubes were set in 
order to analyze the intensity of the emitted fluorescence in five 
non-overlapping bands covering most of the visible light spec- 
trum. Bands were selected to include the acknowledged emission 
peak for the most classical fluorophores (Figure 1). The fluo- 
rescent signal arising from the different fluorophores was then 
collected on these 5 NDDs under two different excitation con- 
ditions, sequentially using excitation wavelengths at 800 nm and 
940 nm. The contribution of each fluorescently labeled compo- 
nent on the 10 NDDs (5 from the 800 nm excitation and 5 
from the 940 nm excitation) is depicted in Figure 2. It should be 
noted that the background signal has to be measured on each 
NDD and subtracted from the fluorophore signal acquired. The 
background signal arises both from electronic and shot noise 
but also from brain autofluorescence. Contribution of the aut- 
ofluorescence to the fluorescence signal depends both on the 
expression level of the fluorescence probe and on its relative 
brightness. It is expected to be larger for thin objects whose 
size is below the voxel resolution. Comparison between raw 
and background subtracted images indicated that endofluores- 
cence accounted for 10-15% of the reference fluorescence signal 
for the dendritic processes and 7-10% for somas in our mice 
(Supplementary Figure 1). The spectral signature of the back- 
ground and of the fluorophores is moreover expected to vary 
according to the imaging depth due to differential scattering and 
absorption. Variations due to depth were limited to less than 15% 
over the first 300 [im below the dura (Supplementary Figure 2). 
Imaging at greater depths would require measurement of the 
spectral signature of used fluorophores at various depths to opti- 
mize spectral unmixing. For each biological component, the 
NDD recording the highest relative signal was defined as the 



Reference NDD for this component: NDD2-800 for vasculature 
(Cascade-Blue), NDD3-800 for astrocytes, NDDl-940 for tumor, 
NDD3-940 for neurons and NDD4-940 for CDllc-positive cells. 
For tumor cells, although the contribution was slightly higher in 
NDD5-940 than in NDDl-940 (33.27 vs. 28.50%, Figure 2F), the 
latter was preferred as it minimized the contaminations by other 
fluorophores. A sixth signal corresponding to the second har- 
monic generation reflected by the dura-mater and peritumoral 
coUagen was analyzed on detector NDD2-940. 

Spectral unmixing was further improved by subtracting the 
signals of contaminating fluorophores on the Reference NDD of 
each biological component (see Methods). For example, the sig- 
nificant contribution of astrocyte signals in the reference NDD 
for vasculature (NDD2-800) was removed (arrows in Figure 3B) 
as well as the contribution of neurons (arrows in Figure 3C) 
in the reference NDD for CDllc-positive cells (NDD4-940). 
Alternatively, deconvolution algorithms available as free ImageJ 
plugins can also be applied (Joachim Walter's Spectral Unmixing 
plugin based on Zimmermann et al., 2002; Neher and Neher, 
2004). 

The combination of multichannel/multiexcitation 2PM and 
spectral unmixing methods therefore facilitates detailed micro- 
metric characterizations of the cellular processes of various cell 
types. In the context of our GBM model, the fine processes of 
CDllc-positive immune cells and connexin-positive astrocytes 
could be independently visualized along the blood vessel walls. 
This supports the viewer that the latter actively contribute to the 
blood-brain-barrier formation (Abbott et al, 2006) (Figure 3D1, 
arrowhead and arrow, respectively). Close apposition and interac- 
tion between CD 1 1 c-positive cells and tumor cell bodies were also 
highlighted inside the tumor core (Figure 3D2) consistent with 
the involvement of CD 1 1 c-positive antigen presenting cells as part 
of the immune response that fights against tumor progression 
(Prins et al, 2011; Gabrilovich et al, 2012). 

The possibility to compute orthogonal reconstructions from 
XY images acquired at different depths in the brain of the anes- 
thetized animal provides access to morphological information 
that would be lost with classical postmortem histology tech- 
niques. In our experimental conditions, the vertical branching 
pattern of blood vessels (Figure 4A, right of the picture) and the 
track used by apical dendrites of layer 4-5 pyramidal neurons 
(Figure 4B, arrows) could both be followed reliably. Something 
that is normally difficult to visualize in fixed thin serial sections. 
It should be noted that acquisition parameters were set such that 
the fine processes of Thyl-GFP cells (neurons) could be readily 
visualized at the expense of cell bodies, which were saturated and 
could therefore not be properly unmixed. Nevertheless, cell bod- 
ies of neurons could be qualitatively differentiated from astrocytes 
somata (Supplementary Figure 3). If required, different acqui- 
sition parameters can be chosen that optimize imaging of cell 
bodies rather than fine processes. The progressive infiltration 
of CDllc-positive cells along the dorso-ventral axis can simi- 
larly be highlighted in orthogonal images (Figures 4C,F, arrows). 
Accumulation of CD 1 1 c-positive cells was pre-dominantly found 
inside the tumor while relatively low numbers of these cells were 
found in the surrounding brain parenchyma (Figures 4E,F). The 
in-depth lateral spatial resolution deteriorates faster inside the 
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FIGURE 3 I Five-color intravital two-photon imaging of glioblastoma 
tumor and its dynamic environment. (A) 5-colors tile scan (3x3 
images) of the GBM and its microenvironment; inset: macroscopic image 
of the cranial window with the area covered by the tile scan (large 
white square) and the area covered by a single plane acquisition (small 
dotted square). (B) Image of vasculature taken from NDD2-800 before 
(B1) and after (82) spectral deconvolution; arrows: location of identified 
astrocytes that are removed from the image after spectal deconvolution; 
inset: zoom (100 x 100 (im). (C) Image of CD11 c-positive cells taken 
from NDD4-940 before (CI) and after (C2) spectral deconvolution; 
arrows: location of identified neurons that are removed after spectal 
deconvolution; inset: zoom (100 x 100 M-m). Note that the background 



generated by DsRed positive tumor cells is removed from the image 
(D). Examples of cellular interactions at the level of their fine processes. 
(D1) zoom of a CDllc-positive cell wrapping its processes around a 
blood vessel (arrowhead) and of an astrocyte whose processes contact 
the vessel walls as part of the blood-brain-barrier (arrow); (D2): zoom of 
CDllc-positive cells surrounding tumor cells in the tumor core. (E) Same 
field of view at day 20 (d20) and day 21 (d21) post-implantation; note 
that astrocytes are stable (arrows) while CDllc-positive cells appear or 
disappear (arrowheads). Colors: blue: vasculature, cyan: astrocytes, green: 
neurons, yellow: CDllc-positive cells, red: tumor cells. Scale bars: 
(A): 200 |Am (inset: 500 |Am), (B-E): 100 iim. All images were taken at 
150 |im depth. 



tumor than in surrounding healthy tissues due to the highly scat- 
tering properties of the dense and structurally disorganized tumor 
tissue (Figures 4A-C,E,F). A more powerful femtosecond pulsed 
laser in less diffusible spectral range (Kobat et al, 2009) (few watts 
and above 940 nm) aiming at efficient two-photon excitation of 
red-shifted fluorophores and triphoton excitation of blue shifted 
ones, might be used to overcome this limitation. 

Spectral deconvolution of the signals from NDD2-940 out- 
lines the superficial collagen fibers (Figures 4A-C) and reveals 
their privileged orientation (Figure 4D). SHG signals from col- 
lagen fibers(Cicchi et al., 2013) of the dura-mater could only be 
visualized in the most superficial 50 \im below the glass coverslip 



(Figure 4D), due to the facts that reflected SHG signals are intrin- 
sically smaller than transmitted ones. Any generated blue photons 
are easily scattered and absorbed by biological tissue. 

The 20X objective used for this study gave access to a field of 
view of 400 x 400 microns. Larger fields of view covering sev- 
eral millimeters were acquired by tiling several elementary fields 
(Figure 3A). Such microscopic imaging of macroscopic fields of 
view provides unique pre-clinical information about the struc- 
tural heterogeneities of tumor vessels and intratumoral immune 
cell distributions relatively to healthy brain regions. Because two 
sets of images have to be acquired at two different excitation wave- 
lengths at a 2hz acquisition rate for 512 x 512 pixels images, 
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FIGURE 4 I Distributions of six components assessed by intravital 
two-photon imaging. (A-C) Orthogonal reconstructions obtained from a 
stack acquired from 0 to 300 M-m below the glass coverslip with a 
z-step of 3yim. Each YZ image shows the maximum intensity 
projection over 10 microns along the X axis. (A) Note the vertical 
orientation of major brain vessels; dotted lines: the levels of xy-sections 
are shown in D-F. (B) Apical dendrites of neurons are visible (arrows). 
(C) Astrocyte (arrowhead) and CDllc-positive cell (arrow) also visible in 



F. (D-F) xy-sections taken at 45 (im (D), 145 iim (E), and 205 iim 
(F) below the glass coverslip. (D) Note the SHG signal (magenta) 
corresponding to the collagen fibers at the level of the dura-mater. 
(E) CDllc-positive cells have invaded the tumor but not surrounding 
healthy tissues. (F) Astrocyte (arrowhead) and CDllc-positive cell 
(arrow) also visible in (C). Colors: blue: vasculature, cyan: astrocytes, 
green: neurons, yellow: CDllc-positive cells, magenta: SHG signal 
(dura-mater), red: tumor cells. Scale bars: 100 |xm. 



approximately 40 min were required to acquire large field of view 
recordings as presented in Figure 3A (9 fields of view with 3 
micron z steps over 400 microns depth). The non-invasiveness 
of 2PM allows for repeated imaging of the same areas in order 
to study tumor progression and micro-environmental dynamics 
on a daily basis. Repeated imaging of tumor burden on day 20 
and day 21 post-surgery (Figure 3E), for example, revealed the 
dynamic distribution of immune CDllc-positive cells (arrow- 
heads) while astrocytes, in contrast did not move (arrows). 

Such information is of interest given the strong influence that 
the micro-environment has on GBM development (Charles et al.. 



2011; Hanahan and Weinberg, 2011). Multiparametric dynamic 
imaging of the tumor and its micro-environment provide unique 
data that could benefit our understanding of current models of 
tumor progression and identify new therapeutic targets. Whereas 
dynamic monocolor imaging allowed for quantification of the 
daily tumor progression rate (Figure 5A1), dual color imaging 
of blood vessels and tumor cells further highlighted the influ- 
ence of tumor growth on the peritumoral vessels (Figure 5A2). 
In the tumor margin, the distances between vascular branches 
increased over time along the vertical axis while they shrunk 
along the horizontal axis (Figure 5B). These changes could result 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



February 2014 | Volume 8 | Article 57 | 6 



Ricard and Debarbieux 



Six-color intravital two-photon imaging 



d20 



If 



d21 



It 



d22 



If 



Distance variations between blood vessels over 
2 days (D20 to D22) 











— ( — 


□ Horizontal ■ Vertical 



Movements 





■ in between blood 
vessels (n=10) 



■ between blood 
vessels and 
astrocytes (n=50) 

between blood 
vessels and neurites 
(n=25) 

between D21 and 
D22 

Distribution of the distance between CD11c-«- 
cells and an invariant vascular landmark 





I day20 
iday21 



□ day22 



O-SOijn 50-100 ijm 100-150 (jm 150-200 urn 



FIGURE 5 I Cell dynamics and deformations due to tumor growthi can 
be tracl<ed over time. (A1) Tumor margin are highlighted at day 20 (d20, 
white small dotted line), day 21 (d21, yellow dotted line) and day 22 (d22, 
green line). Tumor expansion rates were measured at 25 ± 2|im/day. (A2) 
Tumor margins at d20, d21, and d22 were overlaid on vasculature images 
at day 20 and day 22. Notice the decrease of the distances between 
blood vessels in the horizontal direction (in the direction of tumor growth, 
horizontal double-arrow) and an increase in the vertical direction 
(perpendicular to the direction of tumor growth, vertical double-arrow). 
Movements between d20 and d22 are highlighted by the mismatch 
between arrowheads and vessel walls (dotted line on d22). (B) 
Quantification of the distance variation (|im) between blood vessels in the 
horizontal direction (in the direction of tumor growth, light gray) and in the 
vertical direction (perpendicular to the direction of tumor growth, dark 
gray). (C) z-projection of 6 slices over ^5\^m in the same field of view at 



day 20, day 21, and day 22 on a xyz registered image. Vascular landmarks 
(blue, arrowhead), astrocytes (cyan, asterisk) and neurites (green dots, 
arrow) appear stable over the observation period. (D) Quantification of the 
movements (between d20 and d21, left panel and between d21 and d22, 
right panel) between blood vessels (blue, n= 10 measurements), of 
astrocytes relative to blood vessels (cyan, n = 50 measurements) and of 
neurites relative to blood vessels (green, n= 25 measurements). (E) Same 
field of view as in C (z-projection of 6 slices over 15|im) at day 20, day 21, 
and day 22 on a xyz registered image showing only the CDllc (yellow), 
blood vessels (blue) and tumor (red) channels for the sake of clarity. Note 
the dynamics of one CDll c-positive cell identified by an arrow. (F) 
Frequency distribution of the distances to peritumoral CDll c-positive cell 
somas relative to a reference vascular landmark (arrowhead in E) at day 
20, day 21, and day 22. Colors: blue: vasculature, cyan: astrocytes, green: 
neurons, yellow: CDll c-positive ceils, red: tumor cells. Scale bars: 50|im. 



either from vascular remodeling as observed in peritumoral areas 
for other tumor models (Ricard et al., 2013), or from the com- 
pression exerted by the expanding tumor onto the tissue. To 
discriminate between the two hypotheses, we here used infor- 
mation from the color channels corresponding to astrocytes and 
neurons. The displacements of identified astrocytes or neurons 



relative to identified vessel branches exhibited amplitudes sim- 
ilar to the movements of the vessels themselves (Figure 5D) as 
expected from a global deformation of all the cellular compo- 
nents. Over 3 days, deformation was restricted to a few microns 
in amplitude, indicating that the tumor barely influenced the 
structure of the parenchyma beyond 50 microns away from its 
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invasive margin (Figure 5C). Inflammatory cells represented the 
most dynamic features of the pathological brain as illustrated by 
their quickly changing densities and the distribution patterns of 
yellow Cdllc(+) cells (Figure 5F) inside stable neuroglial and 
vascular meshworks (Figures 5C,E). 

Since this approach can be adapted to study other cell popu- 
lations or cell interactions with extracellular matrix components 
using relevant strains of transgenic fluorescent animals, we believe 
that our setup, protocol and animal model open new perspectives 
in the field of neuro-oncology. Generally speaking, the multiexci- 
tation multicolor approach described here can be applied to study 
many pathologies including spinal cord lesion (Fenrich et al., 
2012) but also diseases of other organs such as intestines (Ricard 
et al, 2012a; Zhuo et al, 2012), liver (Honda et al, 2013) or lungs 
(Fiole et al., 2012). Importantly, this multiplexed approach should 
prove valuable for pharmaceutical research since the effect of can- 
didate therapeutical molecules can be simultaneously assessed in 
vivo, on several cellular targets. In addition to a reduction in 
the use of experimental animals, the current methods provide a 
powerful tool for studying drug specificity and the optimal time 
windows for treatment. 
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Supplementary Figure 1 | Astrocytes soma and process spectral 
signatures. Normalized (%) contribution of astrocytes soma (top panel) 
and process (middle panel) on the NDDs for excitations at 800 nm 
(NDD1-800-NDD5-800) and 940 nm (NDD1-940-NDD5-940) without (left) 
and with (right) background subtraction. All measurements were 
performed at lOOiJim below the dura-mater. Measurements were realized 
on 5 different ROI and error bars are s.e.m. 

Supplementary Figure 2 | Astrocytes spectral signature at 50, 150, and 
250 M.m below the dura-mater (A) Normalized (%) contribution of 
astrocytes on the NDDs for excitations at 800nm (NDD1-800-NDD5-800) 



and 940 nm (NDD1-940-NDD5-940) at 50, 150, and 250 n,m below the 
dura-mater without (top panel) and with (bottom panel) background 
subtraction. (B) Measurements were realized on 5 different ROI and error 
bars are s.e.m. 

Supplementary Figure 3 | Neuron and astrocytes somas after spectral 
unmixing. Maximum intensity projection of 5 images with a z-space of 
3 (Am acquired at 200 \Lm below the dura-mater. Spectral unmixing was 
performed with the parameters presented in Figure 2. Acquisition 
parameters were set to highlight the fine neuritic processes which hence 
led to a saturation of the neuron soma. Nevertheless cell bodies from 
neuron (arrow) and astrocytes (arrowheads) can be easily discriminated 
after spectral unmixing. Colors: blue: vasculature, cyan; astrocytes, green: 
neurons. Scale bar: 100ii,m. 
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